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Abstract. The cross section for the dd — > He tttt reaction is estimated near threshold in a two-step model 
where a pion created in a first interaction produces a second pion in a subsequent interaction. This approach, 
which describes well the rates of 27r and r; production in the pd — > '^HeTrTr and dd — > *He?7 reactions, leads 
to predictions that are much too low compared to experiment. Alternatives to this and the double-Zi model 
will have to be sought to explain these data. 

PACS. 13.60.Le Meson production - 14.40.Aq pi, K, and eta mesons 



1 Introduction 

Over the last few years there has been increased exper- 
imental interest in double-pion production near thresh- 
old in several hadronic reactions. These include studies in 
pion-proton ^ and proton-proton collisions ;2 , as well as 
in the pd -^ ^Re tttt pilU^ and dd -^ ''He tttt' P7II5| reac- 
tions. For excess energies Q below about 100 MeV one sees 
no sign of the low mass s-wave tttt enhancement, known 
as the ABC effect Pj, and the maxima in the invariant 
mass distributions tend more to be pushed to the highest 
possible values. 

Due in part to an isospin filter effect, the most spec- 
tacular manifestation of the ABC is to be found in the 
case of dd^ ^Reinr for Q ^ 200 -300 MeV [THj. These 
cross section data, as well those representing the deuteron 
analysing powers jll| . can be well understood within a 
model where there are two independent pion productions, 
through the pn — > dir^ reaction, with a final state inter- 
action between the two deuterons to yield the observed 
a-particle |12| . Since the pn — > dir^ amplitudes are domi- 
nated by p— wave production, driven by the A isobar, this 
leads to much structure in the predictions. Although such 
double-Z\ effects are generally observed in the medium en- 
ergy data |1()II12J . there is little evidence of them nearer 
to threshold "5",T, 13 . Furthermore, the cross sections mea- 
sured at low energies are over an order of magnitude higher 
than the predictions of models behaving like the square of 
p-wave production, where the amplitudes must be propor- 
tional to Q. 
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In an alternative approach to the pd — *■ ^He tttt reac- 
tion, the low energy cross sections have been discussed 
in terms of a two-step model, where a pion is produced 
through a.pp —^ d?:'^ reaction on the proton in the deuteron, 
with a further pion being created in a secondary n^n — > 
Ptt'^tt'^ reaction ^S]. There are, of course, other contri- 
butions related to this through isospin invariance. Semi- 
phenomenological models of the 7r~^n — > mr^Tr^ amplitudes 
show strong s-wave production, behaving rather like a 
contact term, plus another contribution involving the de- 
cay chain 7V*(1440) -> Z\(1232)7r -^ Ntttt 14. The s- 
wave term is sufficient, in the two-step model, to lead 
to reasonable agreement with the available data on the 
pd -^ '^He TTTT total cross section. Moreover, combined with 
p-waves required by the decay chain, it reproduces the 
shift of the mass spectrum away from the ABC region to- 
wards that of higher missing masses I3!,T . It is therefore 
reasonable to ask whether a similar approach could not be 
usefully tried for the low energy dd — > ^He tttt reaction. 

The two-step model with an intermediate pd — *■ '^He tt'^ 
step has in fact been applied successfully to the produc- 
tion of 77-mesons in the dd — > *He 77 reaction near thresh- 
old |15j . where it reproduces reasonably well the magni- 
tude of the total cross section [TCItTT] . The approach is here 
extended in section [21 to describe the dd -^ ''He tttt reac- 
tion, using the same ttN — > ttttN amplitudes as those that 
worked for pd — > '^ He tttt. The other element that is crucial 
for the evaluation of this model is the cluster decomposi- 
tion of the a-particle in terms of ^He n/ '^Hp constituents. 
This is discussed in section|31 where we rely on the work of 
the Argonne group 20 . The results presented in section^] 
show that the model is capable of describing the shape of 
the TTTT effective mass distribution, without the oscillatory 
structure predicted by the double- Z\ model ^2]. However, 
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the total cross section estimates fall over an order of mag- 
nitude below the experimental results found at low ener- 
gies [5irnil3i. These data have low statistics and limited 
acceptance, though they will be supplemented by more 
precise results expected soon from CELSIUS |8j. Since 
neither this nor the double-zi model gets even close to 
the observed production rates, alternative approaches are 
necessary. 



2 The reaction model 



■ ^HeTTTT amplitudes, in 
and TT+n — *■ (7r7r)°p, is 



The two-step model for the dd 
terms of those for pd ^ "^Htt 

depicted in Fig. ^ Contributions involving intermediate 
^He and tt^/tt" are all related to the results for this di- 
agram through isospin invariance. Due to the identical 
nature of the incident deuterons, there is a similar set of 
diagrams where the initial production takes place on the 
upper deuteron. 




Fig. 1. Two-step model for the production of tt^tt" and 
tt^tt" pairs through the dd — > *He tttt reaction. There are also 
contributions related to this by isospin in addition to the terms 
arising from the interchange of the two deuterons. 



we neglect the deuteron D-state and the Lorentz boost of 
the wave functions, this can be written as 
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where the particle masses are denoted by m^ . The reduced 
nuclear matrix element is 
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— = (T • ed {-Apa ■ ed' - iBpd • (cd' x cr)} 
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v2 
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where the (e^, ed') are the polarisation vectors of the two 
incident deuterons and the kinematics are defined as in the 
figure. The 5-state momentum-space wave functions for 
the deuteron and the triton-proton configuration of the 
a-particle are denoted by (pd{<l) and f/^Q, (fc) respectively. 
In the forward and backward (cm) directions, only two 
terms are needed to describe the spin structure of the 
dp — > '^He7r° amplitude. Using two-component spinors to 
denote the ^He {uh) and proton (wp), this reads 



'p I 



M{dp -> ^Hctt") = ul [Apd ■ ed + iBpd ■ (cd x cr)] - 

(2.4) 

where Pd and p^ are the momenta of the incident deuteron 
and produced pion respectively. In our normalisation, the 
unpolarised differential cross section and deuteron tensor 
analysing power i20 are given in terms of the two dimen- 
sionless spin amplitudes A and B by 

da 

dr? 



t2 



1 P-Tl ||2 


+2\B\^' 


V2 


i\B\^-\A\'n 

_|y^|2+2|S|2_ 


1 



(2.5) 



The cross section corresponding to such a diagram 
has been evaluated for the dd ^ arj reaction ^15j and ward 
we follow closely the techniques used there. The unpo- 
larised dd -^ ^He TTTT differential cross section is expressed 
in terms of the Lorentz invariant matrix element A4 through 



and these observables have been well measured in collinear 
kinematics at Saturne [Blj . 

For deuteron kinetic energies of interest here, the back- 



9p^ = 180°) values of ^20 are strongly negative, so 
that 1^1 > \B\. In the 0.5 - 0.8 GeV range the results may 
be represented by 
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Here pd and Pa are the initial and final momenta in the 
overall cm system where the total energy is W. The angles 
f2a are also in the total cm system, whereas the tttt relative 
momentum fc^^ and its angles /2^7r are evaluated in the 
dipion rest frame, where the total energy is rriTrTr- 

The matrix element of Fig. ^ involves the integration 
of the pion propagator between the two production ver- 
tices over the two Fermi momenta k and q). If initially 



where the deuteron kinetic energy Td is measured in GeV. 
The spin structure of the n^p — > 71*^71*^71 amplitude is 
unique near threshold: 

M{-K^p -^ 7r*^7r"n) = a{mT^Tr,Q) u'^cr ■ p^Up . (2.7) 

In terms of the amplitude a, the unpolarised differential 
cross section is 
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Here Pp and p„ are respectively the initial and final nu- 
cleon momenta, Wt^m the cm energy in the ttN system, 
and Q = W-^n — Ztott — rriN the excess energy above the 
two-pion threshold. 

The low energy data in different isospin channels are 
well described by the Valencia model J^ and this allows 
one to project out the / = combination required as input 
in equation H2.3|l . The results can be parameterised as: 



1 



647r3 



|a(m^^, Q)|2 = (1.092 - 0.0211O + O.OOOISQ^) 

-(4.18 + 0.0075Q - 0.00098Q2) x 

-(47.65 - 0.935Q + 0.00743Q^) x^ ^b/MeV^ , (2.9) 



where x = {irLj;^ — 2m^)/mT^. 

Due to small recoil corrections, this parameterisation 
should be used at an excess energy of Q', where 

Q' w xvIt^ + {Q — xmT^){l + 2nnT^ / nia) / {1 + 2raT^/rap). 

(2.10) 

Since large Fermi momenta are not required in the 
two-step model, the dp -^ ^Hctt" and ttN -^ ttttN ampli- 
tudes can be taken outside of the integration in eq. (|2.2|l 
with the values pertaining at zero Fermi momenta. Con- 
sidering only the positive energy pion pole, to first order in 
k and q one is left with a difference between the external 
and internal energies of 

AE^E,^t-Eint^AEa + k-W + q-V, (2.11) 

where 
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Here {Et^, Ed, Et, En) are the pion, deuteron, triton, and 
nucleon total energies, evaluated respectively at momenta 



\Pa - hPd^ Pd^ iPd^ and |p„. 



The relativistic relative velocity vectors V and W de- 
pend only upon external kinematic variables: 
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The resulting form factor 
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then involves a one-dimensional integration over wave func- 
tions in configuration space. In terms of this form factor 
the dd -^ ''He tttt differential cross section becomes: 
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where iV„ is the normalisation of the "^He wave function 
and the extra form-factor contribution coming from the 
interchange of the two incident deuterons is indicated. All 
isospin factors have been included, but it must be stressed 
that in eq. H2.16(l A and B refer to the dp -^ ■^Hctt'^ 
and 7r"p — > n'^Tr'^n charge states respectively. Isospin in- 
variance dictates that tt+tt^ production in dd —^ *He tttt 
should be a factor of two larger than tt^tt", but this sim- 
ple rule is significantly modified near threshold by phase- 
space factors arising from the pion mass difference. 

Two further refinements need to be implemented in 
eq. H2.16|l before comparing its predictions with experi- 
ment. Although the final a-particle is slow in the cm sys- 
tem, relativistic corrections cannot be neglected for the 
incident deuterons. These can be included by boosting Vii , 
the longitudinal component of V, i.e by taking as argu- 
ment of the deuteron wave function J15j 



V = iV^,EdV\\/md) 



(2.17) 



Secondly, the effects of the deuteron D-state have to 
be considered and this can be accomplished by introducing 
two form factors: 

SsMV',W,AEo) 

/•OO 

= 27r2/ dte'*'^^°'F*{-tW)<Ps,D{tV'), (2.18) 
Jq 

where 'Ps.oi''') are the deuteron S~ and Z?-state configu- 
ration space wave functions normalised by 



/■OO 
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The S- and Z?-state form factors enter in different com- 
binations for the A and B amplitudes and, after making 
kinematic approximations in respect of the _D-state com- 
bined with the Lorentz boost, one finds 
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where, as in eq. (|2.16|l . it is assumed that contributions 
from form factors resulting from the interchange p^ <^ 
—Pd have been included. 
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3 The '*He wave function 

Over the last few years there has been remarkable progress 
in ab initio calculations of the structure of light nuclei us- 
ing variational Monte Carlo techniques [20; . Starting from 
realistic nucleon-nucleon potentials, it has been possible 
to identify various cluster sub-structures in nuclei as heavy 
as ^Be. The results for the unnormalised *He : ■^Hp overlap 
function in configuration space are shown in Fig. 12 where 
the error bars arise from the sampling procedure. 



• e°V(0 




.•^♦♦v, 



f} 



ANL ■'Hel'^Hp overlap function 
(points) compared to fit 



r(fm) 



Fig. 2. Unnormalised *He:'^Hp overlap function as a func- 
tion of the ^H-p separation distance. For the purposes of pre- 
sentation, this has been multiplied by re°'^, where the charge 
average a = 0.854 fm~^. The results of Ref. |20| have been 
parameterised as in eq. H.'-{.21^ . 

The overlap function has been parameterised by 



■0(r) 
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^a„e-""^ 



(3.21) 



where a = 0.854 fm ^ represents the average for the '^Hp 
and ^Hen configurations. To ensure good behaviour at the 
origin, the final parameter is fixed by ag = — X]n=o'^"' 
while the other values are sequentially 5.1525, —2.8414, 
-45.1886, 110.7401, and -100.3994. The normalisation 
has been chosen such that 



r^ [tl;{r)f dr = N^ . 



(3.22) 



In the spirit of our approach here to pion production, 
where only these cluster contributions are considered, it 
is appropriate to assume that the p '^H and n '^He compo- 
nents saturate the wave function and take Na = 4 rather 
than the reduced spectroscopic factor obtained in ref. 20 . 



4 Results and Conclusions 

In Fig. O we show the prediction of the shape of the 
missing-mass distribution for inclusive two-pion produc- 



tion at an excess energy of Q = 29 MeV with respect to 
the 27r° threshold. Though the general form is in good 
agreement with the experimental data j(i,:7 , the results 
are too low by almost a factor of twenty! The peak of the 
distribution is predicted to be a little to the right of that 
corresponding to pure phase space, which is also shown. 
Such a feature was clearly observed for the pd — > ^HeTTTt 
reaction at low energies ^3, but the limited statistics in 
the dd case prevents us from drawing firm conclusions 
here. Estimates in the double- Z\ model ^^, which agreed 
convincingly with the data in the resonance region, were 
even poorer compared to the near-threshold data. Apart 
from being a similar factor of twenty too low, this model 
also predicted significant structure in the mass distribu- 
tion which is absent from the experimental data. 




270 280 290 300 

M. [MeV/c') 

Fig. 3. Missing-mass distribution for the dd -^ "^HeX reac- 
tion measured at 570 MeV 0. The chain curve corresponds to 
TT^TT^ production within the two-step model whereas the solid 
one represents the sum of this and tv^tt^ production. The pre- 
dictions are normalised to the integrated cross section by mul- 
tiplying by a factor of 17.6. The dotted and broken curves are 
the similar predictions from phase space, again normalised to 
the total rate. 



The discrepancy is similar for the other low energy 
data 13 , though here the acceptance was small and as- 
sumptions had to be made in order to extract a total cross 
section. In Fig. 21 we show the estimates of the total cross 
sections for the production of charged and neutral pions 
within the two-step model. 

The central problem for any model that attempts to 
describe the dd -^ "'He tttt cross section at low energies is 
that the production of isoscalar pion pairs is very similar in 
deuteron-deuteron and proton-deuteron collisions. Thus 
at Q = 29MeV gE], 



atotidd-^ ''He tttt) _ 40 nb _ 2 
c^tot {pd ^ •^HeTTTr) 60 nb 3 



(4.23) 
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Fig. 4. Total cross section for the dd -^ HevrTr reaction. The 
experimental data from Ref. ^ (star) and Ref. [13] (circle) are 
compared to the predictions of the two-step model scaled by a 
factor of 17.6. The chain curve corresponds to tt^tt" production, 
the broken to tt^tt", and the solid to their sum. 



tive near threshold ;.24-. However, the effect will steadily 
diminish with energy and eventually change sign at the 
resonance. Crude estimates indicate that any effects due 
to such final-state interactions are likely to be less than 
50%, even very close to threshold, and so they are very 
unlikely to provide the explanation of the defect. 

Now, although we have normalised the ^He wave func- 
tion as if it consisted purely oi p^H/n^He pairs, in reality 
the "^H in such a nucleus is on average smaller than the 
physical triton. Nevertheless we have taken the amplitudes 
for pd — > "^He tt*' from the measured data. The same crit- 
icism can be levelled at the double-Z\ model, where the 
final deuteron in the pp -^ d7r+ input would really be 
required for a small deuteron. If there were major correc- 
tions due to such effects they would be likely to be present 
at all energies and hence destroy the excellent agreement 
with data achieved at higher energies J^l- Further inspi- 
ration is therefore clearly needed to resolve this dilemma. 

This work has been much influenced by long-standing discus- 
sions with Pia Thorngren, which have been beneficial to both 
sides. One of the authors (CW) is appreciative of the hospi- 
tality shown to him by Uppsala University. Support from the 
EtaNet programme of the EU is gratefully acknowledged. 



On the other hand, the production of the rj meson is 
much weaker in the dd case, with the ratio of the squares 
of the amphtudes being 



\f{pd -^ ^He?7)| 



1 
50' 



(4.24) 



though perhaps this would be increased by a factor of two 
if corrections were made for the effects of the ?7-nucleus 
final-state interaction. Since the low energy pd — > '^He tttt, 
pd ^ '^Hery, and dd — > ^Hery cross sections are all suc- 
cessfully described by the two-step model, a factor of ten 
undershoot in the dd —>■ * He tttt case is not too surpris- 
ing. The crude comparison made here does not take into 
account fully the spin-parity considerations and the pre- 
diction would have been increased by more than a factor 
of two if the sign of the deuteron _D-wave were reversed 
in eq. ifT^ . 

Given that neither the two-step nor the double-Z\ 
model seems capable of describing the magnitude of the 
dd -^ ^He TTTT cross section near threshold, one must seek 
alternative explanations or modifications to the existing 
mechanisms. Other diagrams, such as that of the impulse 
approximation where the process is driven by pd ^ "^He tttt 
with a spectator nucleon, give very small cross sections due 
to the large momentum transfer. We have not included any 
specific TTTT final-state interaction, but the s-wave scatter- 
ing lengths are relatively small "S^" and, in any case, the 
effects are implicitly included through the use of empirical 
TriV -^ ttttN amplitudes (14) . 

The interaction of the low energy pions with the final 
"^He nucleus might enhance the cross section since it is 
known that the p-wave pion-nucleus interaction is attrac- 
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